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1 Introduction 

Recent years have seen a number of holographic constructions aimed at modeling 
fermions at finite density in 2+1 dimensions (see [1-3] for recent reviews). A primary 
motivation has been to find a suitable gravity dual for non-Fermi liquids, that is 
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metals whose low energy dynamics differ from those of the usual Landau- Fermi liquid. 
In conventional metals the robustness of Fermi liquid theory has been argued in terms 
of an IR stable Gaussian fixed point of the renormalization group governing the low- 
energy excitations of charged fermions about a Fermi surface in 3+1 dimensions 
[4, 5]. If, on the other hand, the physics of the conduction electrons is essentially 
confined to two spatial dimensions this no longer applies. The low energy excitations 
of charged fermions about a Fermi surface in 2+1 dimensions are strongly coupled 
and perturbative field theory cannot be relied upon [6-8]. 

In this context, gauge/gravity duahty may provide useful new insights and con- 
siderable effort has been put into developing holographic duals of strongly coupled 
fermion dynamics in 2+1 dimensions [9-27]. An important step is to demonstrate 
how a Fermi surface may be encoded in a holographic geometry. One approach is to 
couple a probe fermion to the bulk gravitational sector and compute the resulting 
spectral density [9-12, 21]. This has led to many interesting results, including sig- 
nals of Fermi surfaces, Fermi liquids, and even non- Fermi liquid behavior. Interesting 
characteristics of the probe fermion arise because the spacetime geometry and other 
bulk fields serve to provide a strongly coupled gapless sector to which the probe is 
weakly coupled [15]. The probe fermion results are, however, rather indirect from 
the point of view of the holographic background which, by definition, does not know 
about the presence of probes. 

In the present paper, we instead look for signals of Fermi surfaces in static 
susceptibilities, such as a current-current correlator at finite momentum and zero 
frequency [28]. Such correlation functions are inherent to the holographic setup and 
the relevant gauge couplings are determined by the same bulk Lagrangian as the 
gravitational background itself. Specifically, we expect to see non-analytic features, 
so called Friedel oscillations, at twice the Fermi momentum in zero temperature static 
current-current correlation functions when there are a discrete set of Fermi surfaces. 
Friedel oscillations are expected whenever there is a sharp Fermi surface, and do not 
rely upon the assumptions of Fermi liquid theory. 

The holographic models that we consider have three key ingredients in common. 
One is that the gravitational dual spacetime is asymptotically AdS, which corre- 
sponds to an underlying conformal symmetry in the boundary field theory. Relaxing 
the condition of conformal symmetry is of considerable interest but this will not be 
pursued here. Another shared ingredient is a U(l) electric flux threading the asymp- 
totic AdS region, which is the holographic manifestation of a non-vanishing chemical 
potential and charge density in the boundary theory. The finite charge densities at 
zero temperature indicate that the dual field theory also possesses Fermi surfaces by 
a generalized Luttinger's theorem [29]. Finally, we have chosen to focus on models 
where the asymptotic electric flux is sourced by a finite density of charged fermions 
inside the bulk rather than emanating from behind the event horizon of a charged 
black hole. 
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One such system is an electron star, first considered in [14, 17]. There the ap- 
proach is to increase the charge density of bulk fermions, such that they may be 
considered as a fluid of degenerate fermions. In this limit, there is a continuum of 
bulk Fermi surfaces, distributed radially. Further, in a certain parameter regime, 
one can ignore the effect of gravity on the fluid's local equation of state. By setting 
up the corresponding Oppenheimer-Volkov equations it is possible to numerically 
determine the spacetime's geometry. Electron stars are non-singular geometries pos- 
sessing Lifshitz scaling regimes in the deep interior. One of the initial hopes for 
electron stars was that, by working with macroscopic bulk charge densities, they 
would make fermionic properties visible at the classical level. Indeed, in [21, 23], 
it was shown, at the probe fermion level, that electron stars support many discrete 
Fermi surfaces which satisfy a Luttinger count. On the other hand, as discussed 
in [19], most of the bulk Fermi surfaces do not contribute to de Haas- Van Alphen 
oscillations in a background magnetic field. These are dominated by a single Fermi 
surface while the contributions from the remaining local Fermi surfaces in the bulk 
fluid wash out. We will find similar smearing when we look for Friedel oscillations in 
electron star backgrounds except in this case the signal is completely smeared out. 

Another model for holographic Fermi surfaces, which we will refer to as an "AdS 
hard wall," was introduced in [25]. This model was intended as a toy model where 
a finite number of discrete Fermi surfaces could be made manifest. Interpreting the 
difficulties encountered with extremal black holes and electron stars as a signal of too 
many IR degrees of freedom, [25] instead works within a confining dual gauge theory, 
achieved via a truncation of an AdS geometry by introducing "hard wall" boundary 
conditions at a finite value of the radial variable. The profile of a bulk gauge field in 
the presence of a filled bulk Fermi surface is solved for numerically, neglecting any 
gravitational back reaction, and then the effect of the resulting gauge field on the bulk 
fermions is solved for. These steps are iterated towards a self-consistent configuration 
of the fermions and gauge field. A Luttinger count indicated that the resulting 
theory could describe either Fermi or non- Fermi liquids depending on the gauge field's 
boundary conditions at the hard wall [25]. Morally, this configuration is similar to 
that of an electron star with a non-vanishing charge density carried by bulk fermions, 
but differs crucially in the number of IR degrees of freedom. Also, the character of 
bulk Fermi surfaces is very different than in the electron star geometry. In the 
hard wall construction, a bulk Fermi surface is not a local concept, but instead each 
Fermi surface is assigned a radial mode number. The non-locality in the bulk Fermi 
surface leads to a certain factorization in the way radial and boundary information 
is encoded, implying that the AdS space outside the hard wall acts more or less like 
a finite size box, as discussed in [25]. 

Our understanding of the application of holography to fermionic systems at fi- 
nite density is in its infancy. The electron star and hard wall geometries provide 
two systems through which we may develop our intuition. By contrasting the two 
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backgrounds we compare the effects of bulk Fermi surfaces wliicli are local vs. non- 
local and continuous vs. discrete. In order to see signatures of Fermi surfaces in 
current correlation functions, which involve local bosonic operators, it is necessary 
to include the effects of bulk fermionic loops. This is because fermionic operators do 
not have expectation values. The primary result of this paper will be to demonstrate 
that bulk Fermi surfaces in an AdS hard wall geometry can induce boundary Friedel 
oscillations, while a similar calculation for electron stars shows that the oscillations 
are washed out due to the continuum of bulk Fermi surfaces. These results suggest 
that discrete bulk Fermi surfaces which are non-local in the radial direction induce 
sharp, discrete boundary Fermi surfaces. 

The paper is outlined as follows. We begin with a review of Friedel oscillations 
in a 2 + 1 dimensional perturbative system in Section 2. In Section 3 we briefly 
discuss the AdS hard wall background introduced in [25] and then we identify the 
source of Friedel oscillations. In Section 4 we set up an analogous calculation in an 
electron star geometry but find that there are no boundary Friedel oscillations due 
to the continuum of bulk Fermi surfaces. Finally, we conclude with a discussion in 
Section 5. In order to streamline the presentation the technical details are relegated 
to appendices when possible. 

2 Friedel Oscillations 

In systems with sharp Fermi surfaces, zero temperature, static response functions dis- 
play spatial oscillatory behavior, called Friedel oscillations.^ Heuristically, this may 
be understood as follows: Having a sharp, non-analytic, cut in the quasi-particle 
density of states leads to non-analytic features in momentum space static response 
functions. Upon Fourier transformation, these non-analyticities give rise to oscilla- 
tory spatial features at wavenumber 2kf.'^ Such oscillations occur in all known models 
with sharp Fermi surfaces, for essentially any static response function considered {i.e. 
density-density, spin-spin, etc.). Heuristically, the character of response functions 
changes depending on whether the exchanged momentum "fits" in the Fermi sphere. 
In this section we will outline how this works for 2-1-1 dimensional non-relativistic 
fermions, leaving the details for Appendix A. Understanding the perturbative calcu- 
lation will be useful to help set notation and because it has close parallels with the 

^Strictly speaking Friedel oscillations refer to spatial oscillations occurring in charge density 
perturbations. In this article we will use the term more generally to refer to spatial oscillations in 
generic transport coefficients due to a sharp Fermi surface. 

simple example is given by the Fourier transform of a one-dimensional sharp bump function 
in momentum space of width w and unit height, 

F.T.[B{k)] =sm{wx)/w. (2.1) 
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holographic calculation. Similar effects are expected in non-Fermi liquids since it is 
the sharpness of the Fermi surface that matters rather than the detailed dynamics. 



2.1 Perturbative Friedel Oscillations in 2 + 1 Dimensions 

Friedel oscillations are easily visible in the density-density polarization {i.e. suscep- 
tibility), denoted by 11. For simplicity we will focus on a system of spinless fermions 
with a Fermi momentum, kf. For weakly coupled fermions this is textbook material, 
see for example [30, 31]. 

In a unitary system we can use a Lehmann representation to capture the full 
frequency dependence of two point functions. The time ordered Green's function for 
free fermions in a translationally and rotationally invariant system takes the form 

X 0{\k\-kf) e(kf-\k\) , , 

in (cu, k) space. Using the fermion's Green's function, an application of Kubo's 
formula indicates that the leading contribution to the polarization is 

n(g, z/) = -le^ j j^^^g{k, uj)g{k + q,uj + p). (2.3) 

In the static limit the polarization tensor is real 

lmn(g,0) = 0, (2.4) 

Ren(g,0) = -^V [ dkdO e{kf-k)—^^—. 

2vr^ J hk+q — tjk 



(2.5) 



For spherical Fermi surfaces the integrand is singular when q = 2k f. Even without 
knowing the precise form of the dispersion relation, the polarization tensor is non- 
analytic in momentum space and thus has oscillatory behavior in position space. 

For concreteness, consider non-relativistic free fermions with Ek = In this 
case (2.5) reduces to 

m«) = ~vjikmeik,-k) ^^j^l^ (2.6) 




(2.7) 



The Kubo formula relates the charge density susceptibility to the polarization 
tensor 

x{q) = {pi-q)p{q))R = -^niq) = -n(g) . (2.8) 

In the last equality we have used the fact that the static polarization tensor is real 
to equate its time-ordered and retarded incarnations. In momentum space the non- 
analyticity at g = 2A;j is clearly visible in Figure 1. 



- 5 - 



m 




Figure 1. The non-relativistic static susceptibility, x{q) = ~n(g), in 2 + 1 dimensions for 
spinless fermions of mass m and charge e. 



3 Friedel Oscillations in an AdS Hard Wall Geometry 

As mentioned in Section 1, a number of holographic models have been introduced to 
describe Fermi surfaces at strong coupling. To reduce the number of IR degrees of 
freedom, the author of [25] proposed to work with confining boundary field theories. 
As is familiar from using holography to study QCD, we can model confinement with 
a hard wall placed in the bulk [32]. Although this model does not include the effects 
of a bulk Fermi surface on the metric degrees of freedom, the Fermi surface does back 
react onto the gauge field's profile. Therefore, in this model one can study the effects 
of a bulk Fermi surface outside of a probe limit on fermionic or current correlators, 
but not on stress tensor correlators. Here we will focus on current correlators and 
show that they display boundary Friedel oscillations. 

3.1 Preliminaries 

Following [25], we begin with pure AdS4, 

ds^ = — {-dt^ + dz^ + dx^) . (3.1) 

In order to control the number of fields in the dual field theory's IR limit, one 
introduces a hard wall cutoff, Zm- We can imagine that the cutoff is an approximation 
to some more interesting geometric effect [33] . On this background geometry we place 
both fermions and gauge fields governed by the standard action 

+ W (0 + m) ^ 



4e2" 



S = d X\f^^ 



(3.2) 



The covariant derivative is = + \^ahiiX"'^ — il^fi ■ We will be interested in 
constructing the background gauge field profiles and understanding hnear response 
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theory in the gauge sector on large length scales. For these tasks it is sufficient to 
integrate out the fermions and use the resulting low energy effective action 

S = J rf^Xv^^F^ + Trln + . (3.3) 

In addition to the differential equations coming from Eq. (3.2), one must also specify 
boundary conditions at the hard wall. Because the geometry is regular, it is natural 
to impose self-adjoint boundary conditions on the fermions [25]. Such boundary 
conditions are quite strong, however, ruling out any dissipation in the fermionic 
sector in the large N limit. 

With the boundary conditions in place, we seek to ffil fermion states up to the 
ground state energy as well as determining the self-consistent gauge field proffie. 
In practice this may be done through an iterative process. As discussed in [25], 
Neumann boundary conditions for Aq at the hard wall, are a necessary condition 
to have a boundary Fermi liquid in this set up. Any other self-adjoint boundary 
condition leads to non-Fermi liquid behavior. The resulting background consists 
of a ffiled sea of fermion single particle wave functions, each indexed by boundary 
momentum as well as a radial mode number. In addition one must include the 
backreaction of the Fermi sea on the gauge field. For chemical potentials above the 
gap, the system has compressible [34] gapless excitations and a sharp Fermi surface. 

It is useful to consider the time ordered Green's function of these normalizable 
fermions. Using a Lehmann representation we have 

where the wave functions and energies are those numerically determined when solving 
for the background, and Ei{k) is the eigenvalue of the Dirac Hamiltonian with spatial 
momentum k and radial mode number i. 

3.2 Mean Field Susceptibility 

Here we turn to the problem of linear response theory in hard wall AdS in the radial 
gauge, Az = 0. The (static) classical wave equation for the potential is 

{dl-kl)Ao{k,,z) = 0. (3.5) 

For definiteness we have taken the momentum to be in the x direction and we focus 
on a gauge invariant sector of perturbations AQ{kx). Since the hard wall is a regular 
point in the geometry, it is natural to impose the same boundary conditions that 
are used in constructing the background. We can find the classical bulk to boundary 
Green's function [i.e. the solution which goes to one at the boundary) analytically. 

Go^{kx, z) = cosh.{kxz) — tanh{kxZm) smh.{kxz) , (3.6) 
{p{-kx)p{kx))o = -kxtanh{kxZm) (3.7) 
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As expected, the classical contribution to the susceptibility does not know about 
the Fermi surface. We must go away from mean field theory and consider virtual 
particle-hole exchanges. 

3.3 Polarization Corrections 

In order to proceed, let us imagine integrating out the fermions in the bulk. Graph- 
ically, the first correction to the gauge field's effective action comes from a bubble 
diagram with fermions in the loop.'^ If we expand out the fermion's functional de- 
terminant to second order in gauge fields, we have a term of the form 

Sp^i = — dzdz'dx A^{z, -k)W{z, z', k)A^{z', k) . (3.8) 

The one loop vacuum polarization tensor is written 

/7 72 
^^Tr [M^G{uj, p, z, z')M''G{uj + u,k + p, z' , z)] , (3.9) 

where = — gF^F^. The reason that ^ —qV^ is that we have chosen to work 
with the Green's function involving Xa:/(-2)xI ti^') rather than Xk,i{z)Xk,i{z')- what 
follows we will limit our discussion to the static polarization tensor, n^'^(z, z' , 0,p) = 
U^^lz, z', p). In the background geometry there are three parameters, namely /i , g , e, 
and only two of them are fixed by the classical geometry [25]. The relative size of 
the loop correction (3.8) is determined by the third one.^ 

In general, evaluating the vacuum polarization in curved space is a rather difficult 
problem. There are the usual ultraviolet divergences occurring due to large i and 
k. Also, one might worry about the cornucopia of fields predicted by string theory 
predict in an eventual top down approach. As discussed in [30], we can split the 
vacuum polarization into two pieces: the (uninteresting) contributions present in 
vacuum and the relative polarization due to a chemical potential, 

W'" = n>'''Uac + nz. (3.10) 

If the chemical potential is not too high, we can assume that the vacuum part only 
contributes to the ultraviolet renormalization of parameters in Eq. (3.2), and is not 
relevant to the long wavelength physics associated with Friedel oscillations. 

^At zero momentum the effective action for Aq also has a tadpole contribution analogous to the 
one used to determine the gauge field's background profile. This term does not affect our results 
below. 

^To anticipate the notation used in Section 4 below, we have rescaled the gauge field A^, ^A^, 
as in Eq. (4.2), with L the curvature radius, and k the Newton constant. With this normalization 
and with the metric as in Eq. (3.f), the background action has an overall factor of which we 
have factored out in (3.8). 
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If we denote the vacuum time-ordered bulk fermion Green's function by G^, the 
relative polarization tensor may be written as, 



■p 



Tr M>'G{co,p,z,z')M''G{LO,k+p,z',z) 



(27r)3 



Tr M^G" 



{uj,p,z,z')M''G%u,k + p,z\z) . (3.11) 



We analyze the relative polarization in Appendix B. The full result is in general quite 
involved but fortunately we do not need it in full detail. Let us list some important 
features. First of all, the relative polarization is finite in the bulk UV. This is 
analogous to the usual finite temperature field theory statement that the only UV 
divergences are vacuum divergences. Second, the integrand cleanly splits into terms 
unrelated to the Fermi surface which are analytic in the external momenta near 2kf 
and contributions that have singular points at k = 2k f. These singularities do not 
lead to actual divergences in the polarization but they are the source of bulk Friedel 
oscillations in much the same way as in the non-relativistic example in Section 2.1. 

3.3.1 Illustrative Limits 

As discussed in detail in Appendix B.l, there is a limit where the relative polarization 
radically simplifies and one can obtain an (almost) analytic expression. If we consider 
a high confinement scale (small Zm) then all contributions to the loop diagram with 
non-zero relative radial momentum^ are suppressed by factors of Zm and may be 
dropped. The integrand then only receives contributions from near a Fermi surface. 

We can make a stronger statement if we take a non-relativistic limit and tune the 
chemical potential towards the fermion's bulk gap energy, leading to a Fermi surface 
of small volume. For simplicity, consider a background with a single Fermi surface, of 
mode number 1. As seen in Figure 2, the numerically determined dispersion relation 
is very well approximated by a quadratic form near the Fermi surface 



Assuming this functional form, the vacuum polarization is well approximated by 



as discussed in Appendix B.l. We have introduced an effective expansion param- 
eter, A = which encodes the shape of the spectrum. The source of bulk (and 

^More precisely contributions where the fermion propagators have different radial mode numbers. 



Ei{k) ~ ai + bik'^ . 



(3.12) 




(3.13) 
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0.01 



0.02 




k 



Figure 2. The energy dispersion relation for the i = 1 family of eigenvalues for 
{m,q,e, jd, Zm) = (1, 1, \/3, 6.287, .5) with kj ~ .20 as well as the best fit curve Ei ~ 
-.003 + .8/^2. 

boundary) Friedel oscillation is now explicit in this non-relativistic limit. Interest- 
ingly, there is an important difference compared to the perturbative calculation in 
Section 2.1 in that the polarization is a non-diagonal tensor. This amounts to non- 
vanishing mixed density- current correlation functions in the holographic dual. 

3.4 One-Loop Corrected Bulk-Boundary Green Functions 

Returning to the gauge invariant sector, it is straightforward to set up of the diffeo- 
integral equations for the loop corrected gauge field equations of motion. 



{dl - kl) A^{K, z) = -e^ / dz' [noo (^, z\ K)Ao{K, z') + Yill{z, z' , K)AAK. z')] . 



The absence of metric factors on the left-hand side is a peculiarity of the gauge field 
equation of motion in AdS4. We will work in a background where the gauge non- 
trivial part of Ar^ vanishes such that we can drop the contribution from the second 
term on the right hand side of (3.14). 

If we have a high confinement scale and a small Fermi surface volume, the diffeo- 
integral equation can be solved perturbatively in A for the bulk to boundary Green's 
functions 



C^'^K, z) = G^'^{K, z) + XGf'\K, z) + .... {t - component) . (3.15) 



This expansion should be valid as long as loop effects are small. Expanding (3.14) 
we find. 



(3.14) 
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Figure 3. Correlation functions for {m,q,e, /j., Zm) = (1, 1, \/3, 6.287, .5) with kf ~ .20 in 
a hard wall geometry. On the left, the full density-density correlation function computed 
by including the one- loop effect induced by the polarization term (3.8). On the right, only 
the one- loop contribution to the density-density correlator is plotted. 



The leading contribution, Gq , was determined in Eq. (3.6). Since the leading 
contribution satisfies the physical boundary condition, G^^{k, 0) = 1, the corrections 
must satisfy Dirichlet conditions near 2; ~ 0. 

3.4.1 Numerical Results 

The perturbative equations (3.17) can be solved numerically and using the AdS/CFT 
dictionary it is then straightforward to translate the bulk to boundary Green's func- 
tions into correlation functions in the boundary theory. Figure 3 shows the resulting 
density-density susceptibility with and without the uninteresting classical contribu- 
tion. At the 0{1/N) level we clearly see a non-analyticity at A; = 2kf, the hallmark 
of a sharp Fermi surface. The analysis in this section is performed for Neumann hard 
wall conditions. It would be interesting to extend the analysis to other boundary 
conditions. Intuition from field theory suggests that Friedel oscillations should still 
be visible whenever there is a sharp Fermi surface regardless of the stability of the 
quasi-particles at the Fermi surface. 

The loop corrections involve the effective loop parameter, ^, encoding 1/A^ cor- 
rections which did not play a role in when determining the classical gauge field 
background. The value of the loop parameter does not affect the position of the non- 
analyticity in (3.14) and we have set it to one in order to maximize the amplitude of 
boundary Friedel oscillations in our numerical results. 

This example illustrates several important points relevant when trying to con- 
struct boundary Fermi surfaces from Fermi surfaces in the bulk. Firstly, even when 



0={d!-kl) G^'{kx,z), 
= (9f - kl) Gf^ikx, z) - 




(3.16) 



(3.17) 
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the role of the bulk Fermi surface is manifest in the construction of the background 
profiles, classically computed linear response susceptibilities typically do not reveal 
the surface. As in perturbative descriptions, one is looking for evidence of particle 
hole pairs in bosonic observables, which are only present in loops (though the loops 
are now in the bulk). As discussed in [25], Luttinger's theorem fixes the boundary 
value of kf to be the same as the bulk value. The amplitude of Friedel oscillations is 
thus dependent, but not the wave vector at which they occur. 

Secondly, the presence of a bulk Fermi surface does not guarantee a clear signal in 
the boundary theory. To be visible, effects of bulk Fermi surfaces must act coherently 
across the radial direction. For example, when we compute corrected Schwinger- 
Dyson equations, there are convolution integrals in the radial direction. In order 
to prevent the convolution integrals from washing out the effects of a bulk Fermi 
surface, the spectrum of radial momenta must be gapped. In the model of [25], it is 
the confinement scale which introduces the gap. 

Thirdly, even if one can neglect smearing due to convolutions, it is necessary 
to have factorization between the dependence on the boundary momenta and the 
radial positions. In this hard wall example this feature is visible in the form of Eq. 
(3.13) (or more generally in Eq. (B.3)). The relevant singularities in the integrand 
are only functions of boundary momenta. In the hard wall example it is the radial 
non-locality of the Fermi surface which makes this happen. 

In the next section we will discuss electron star geometries which also have a 
finite density of fermions in the bulk. However, this geometry differs in two crucial 
ways from the hard wall model. The spectrum of radial momenta is not gapped and 
the bulk fermionic features (such as the Fermi momenta) are local functions of the 
radial positions. 

4 Absence of Friedel Oscillations in Electron Stars 

Another holographic model for the physics of finite density fermions are electron 
stars [14, 17]. The basic elements include a fluid of charged electrons treated as 
an ideal fluid of non-interacting particles at zero temperature, i.e. a Thomas- Fermi 
approximation. The star geometry carries charge more efficiently than an extremal 
black hole and hence it is thermodynamically preferred at zero temperature. In 
the deep interior the gravitational solution goes over to a Lifshitz scaling solution. 
In [26, 27], the finite temperature extension of this geometry was found. As for 
holographic superconductors, at low temperatures it is thermodynamically preferred 
for the charged black holes to expel charged matter. The preferred low temperature 
geometry consists of a charged black hole surrounded by a fermion fluid of the same 
sign charge. Above a critical temperature there is a third order phase transition from 
an electron cloud over to an AdS-RN black hole. We will refer to the zero temperature 
geometry as an electron star and to its finite temperature generalization as an electron 
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cloud. Although both the hard wall fermion background and the electron star have 
bulk Fermi surfaces there are important differences, which will allow us to test the 
lessons mentioned in Section 3.4.1. Ultimately, we find that bulk Friedel oscillations 
will wash out, leaving no sharp signature in the boundary current- current correlation 
function in the electron star. 

4.1 Electron Star Background 

The electron star geometry is given by^ 

7-2 

ds^ = — (-f(z)df + g(z)dz^ + dx^) (4.1) 
A = —h{z)dt, (4.2) 

where e is the fermion charge, L is the curvature radius, and k, is Newton's constant. 
In each local patch there is a degenerate fluid of free fermions modeled as an ideal fluid 
described with velocity, pressure, energy density, and charge density, {u^,p, p,a)J 
In the scaling regime the geometry takes the form 

1 hoo 
f ~ z2{s-i) ' 9 ^ 9oo, and ^ ~ ^ • (^•4) 

In the rest of the spacetime the background profiles are non-singular functions of 
the radial coordinate interpolating between an interior Lifshitz region and a near- 
boundary AdS region. 

One of the initial motivations for the electron star model was to make the 
fermionic features more dominant by working at finite fermion densities. Indeed, 
in [19], it was shown that the thermodynamic potential displays de Haas- van Alphen 
oscillations due to a bulk Fermi surface. However it is also clear that the fermionic 
character of the electron star geometry is quite different than in the hard wall model of 
Section 3. Most importantly, there is a continuum of Fermi surfaces in the Thomas- 
Fermi approximation. At each radial position there is an effective local chemical 
potential 

zh 

A=^. (4.5) 

^Note that we have shghtly redefined the metric functions / and g relative to [17]. We will also 
denote the dynamical critical exponent by s rather than the usual z to avoid confusion with the 
radial coordinate. 

''The fluid variables are rescaled as in [17], 

Furthermore, the mass of the fermions and the constant proportional to the spin in the fermion 

2 

equation of state are rescaled according to m = ^rh, /3 — -^^j^P- 
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A second important difference is the electron star's lack of confinement and a relative 
surplus of infrared degrees of freedom. The emergence of a Lifshitz scaling region in 
the geometry's interior implies that there is no gap in the allowed radial momenta. 

In the rest of this section we discuss Friedel oscillations in the electron cloud 
geometry. We work with the finite temperature electron cloud geometry rather than 
the zero temperature electron star. This is purely for convenience when carrying 
out numerical calculations. Since the electron star geometry, including the Lifshitz 
scaling region, is recovered when the temperature of an electron cloud geometry is 
taken to be very low compared to the scale set by the chemical potential, we can 
expect to recover electron star correlation functions from the ones obtained from 
electron clouds in the low temperature limit. Furthermore, as the temperature of 
an electron cloud geometry is given by the Hawking temperature of the black hole 
horizon below the electron cloud while treating the electron cloud itself as a zero 
temperature fluid [26], any bulk "medium" effects should be visible even at finite 
temperature. 

As a second technical simplification we will not look for Friedel oscillations in 
the sound channel, the natural generalization of Section 3.2. Instead, we will use the 
fact that any channel is expected to display Friedel oscillations in the presence of a 
sharp Fermi surface and work with the much simpler shear channel. 

4.2 Electron Star Static Correlation Functions: Shear Channel 

Not surprisingly, analyzing the linear response theory for electron stars and electron 
clouds is considerably more complicated than for the hard wall. If we consider 
a general perturbation with boundary frequency and momentum {co,ky), we can 
classify the allowed perturbations according to their behavior under the residual 
Z2 symmetry, (x — )■ —x) as discussed in [35, 36]. In asymptotically AdS4, there 
are two families of perturbations in the presence of the fluid in the radial gauge 
6A, = 0, 6g,^, = 0, 

odd : {6g^y, 6g^t, SA^, 6u^) 

even : {Sg^^:^, 6gyy, 6gu, 6gty, 6Ay, SAt, Sut, Su^, 5uy). (4.6) 

The family of odd perturbations give rise to shear modes while the even family 
corresponds to sound modes in the dual field theory. Here we will focus on the shear 
channel because it includes a smaller set of fields. 

4.2.1 Static Limit 

In general, the gravitational response theory is complicated by the fact that there are 
gravitational gauge symmetries in the form of diffeomorphisms. In the static limit 
there is a simplification and we only need to retain 5(7*^, 5uxi and SA^. It is straight- 
forward to set up the classical fluctuation equations describing gauge invariant shear 
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modes in the uj limit, 



= {6g 



2z 

h' 



w 

h' 
7 



= {6A^r - gkl6A^ + '-^{Sg'J - ( ^ + ^ ) + '-^6u^ . 



Z2 



(4i 



In addition to the equations of motion, the divergence of the stress-energy tensor 
fixes 



5Ar^ = —fi6u^ 



(4.9) 



For an electron cloud geometry, it is straightforward to compute Euclidean boundary 
susceptibilities by imposing regularity conditions at the horizon. As in the non- 
relativistic example we find a real static susceptibility and hence the Euclidean result 
equals its Lorentzian counterpart. 



2.0 





E. C. {T/T,=.03) 




B. H. iT/T,=.03) 




E. C.{T/T,=.5l) 




B. H. (T/Tr=.51) 




Figure 4. Classical contribution to static correlation function in the presence of a fluid with 

55 ^ 

100' 



parameters (/3 = 10,m = ^) and for RN black holes at moderate and low temperatures 



in units of T^. 



As in the AdS hard wall example we need to consider vacuum polarization effects. 
Indeed, neither the differential equation (4.7) nor the boundary conditions can know 
about a Fermi surface. Within classical fluid response, there is nowhere for non- 
analyticity in the boundary momenta to enter the problem, and, accordingly, there 
is no evidence of a sharp Fermi surface in Figure 4. A similar result was found in [28] 
for the DA - D8 - D8 system. 



- 15 - 



4.3 Polarized Electron Star 

We would like to include polarization effects on the shear modes' equations of motion. 
This represents a deviation from the idealized perfect fluid modeling. Even in flat 
space it is a complicated problem to go from a microscopic field theory description to 
a long wavelength hydrodynamical picture. In the case at hand we have the added 
challenge that the radial and boundary momenta enter the problem at different levels. 
Boundary momenta are essentially arbitrary, being fixed by boundary conditions, 
whereas bulk momenta are set by the scale of variation in the bulk space time. Here 
we will only be interested in effects which lead to non-analyticity and bulk Friedel 
oscillations, and will ignore other, possibly larger, effects of interactions that only 
contribute analytic terms to the polarization tensor. 

Noting that the vacuum polarization is a short distance effect when compared to 
the scale of variation in the radial direction, we will introduce a poor man's modeling 
of the polarization, 

Spoi = ^ J dzdz'd'k./\g{z) I \g{z') \6A^iz, -k)U^J^ ,,ci^, z', k)5A,{z' , k) (4.10) 

where H^'^ is a radially coarse grained polarization tensor. To obtain an approxi- 
mate value for TI'^q ^ci start from the polarization tensor for relativistic fermions 
in each local Lorentz frame, 11'^'^. Inside the fluid, we will only be interested in the 
behavior of the polarization tensor at long wavelengths in the radial direction. We 
define a coarse grained polarization tensor in flat space by simply projecting on to 
infinite wavelength, 

^'cGJlatik., k^.L.) ^ n^^(0, kL.L.) , (4.11) 

^U^Jaj,Jz,z',kL.L.) = 6{z-z')U^''{0,kL.L:), (4.12) 

where the spatial momentum in each local frame is denoted /cl.l.. 

Locally, we then approximate the (coarse grained) curved space polarization 
tensor by pulling back the (coarse grained) flat space expression, using the local 
values of the chemical potential and Fermi momentum. This last step introduces a 
subtlety by redshifting the spatial momentum which is to be used in the polarization 
tensor. Because the fluid in each patch naturally lives in the local Lorentz frame, 
the momentum flowing through the polarization tensor is that of the local Lorentz 
frame, kL^L = zk. Pulling all this together, we arrive at the polarization contribution 
to the effective action 

Spoi = Y y dz^/\g\^Mz^ -k)W{Q, zk)6A,{z, k) , (4.13) 

written in terms of the flat space polarization tensor evaluated in each fluid patch. 
The full expressions are quite lengthy and may be found in Appendix C. 
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4.3.1 Lack of Friedel Oscillations 



As in Section 2.1, we are only interested in the additional polarization due to the 
finite density of fermionic matter. For brevity we will suppress all contributions not 
relevant for triggering Friedel oscillations. In each local patch, relativistic Friedel 
oscillations are caused by a mildly singular term (see Appendix C for details). 



n'^"(0,fcL.L.) ~ Ar'^'^(A:L.L.)ln 



L.L. 



2ki 



^L.L. 



+ 2k 



+ ..., 



(4.14) 



where N'^'^ has smooth dependence on the momentum near kj and the neglected 
terms are also smooth at kf. With this, the relevant term in the gauge field effective 
action becomes 



S 



Pol 



J dz^/\g\5A^{z,-k) 



zk 



zk + 2 



5A^{z,k). 



(4.15) 



There is an important difference between this polarization tensor and that of Sec- 
tion 3.3. In the hard wall example, the source of the bulk Friedel oscillations was 
z independent. This led to a factorization of radial and boundary information such 
that the bulk Friedel oscillations were coherent across the radial direction. For the 
electron star, the bulk Fermi momentum is a local concept. Even in the deep inte- 
rior of an electron star, where the local chemical potential goes to a constant, the 
polarization is z dependent. 



S 



Pol 



dz^/\g\6A^{z, -k) 



large z 



N'^'izk) In 



zk - 2^ 



2 

oo 



SA^{z, k) . 
(4.16) 



We see that different bulk oscillations are triggered at each radius. 

The absence of boundary Friedel oscillations in electron stars has two basic 
causes: there is a continuum of bulk Fermi surfaces and each surface is localized 
in the bulk. These facts are manifest in the smearing in two ways. First, the argu- 
ment of the polarization is a function of the background profiles for the metric and 
gauge field. Given this, one can expect a degree of smoothing set by the sharpness 
of the transition from the interior Lifshitz to the asymptotic AdS behavior. The 
second, and more important, manifestation is the momentum redshift factor which 
persists all the way into the Lifshitz scaling regime. This is analogous to the smear- 
ing described in [19]. Since the fate of a continuum of bulk Fermi surfaces may be 
of relevance for holographic model building, it is important to understand to what 
extent the two sources of smearing we have identified in this work are independent. 
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One might complain that the polarization effects were included in an ad hoc 
manner in this section. However, the approximations we have made are more or less 
optimal for seeing Friedel oscillations and yet there is still smearing due to the infinite 
number of bulk Fermi surfaces. A more thorough analysis would keep track of the 
radial scale over which the fermions are integrated out, almost certainly leading to 
even less sharp features than seen with the simple projection used here. 

5 Discussion 

In this paper, we have initiated a study of Fermi surfaces in holographic metal- 
lic states, without reference to probe fermions. Specifically, we use static current- 
current correlation functions to characterize boundary Fermi surfaces, focusing on the 
question of when bulk Friedel oscillations trigger boundary Friedel oscillations. We 
discussed two systems of charged bulk matter, not hidden behind horizons: fermions 
in a hard wall truncation of AdS and an electron star geometry. For both of these 
systems a Luttinger count indicates that the dual theory has filled Fermi surfaces 
and for both systems it is necessary to consider loop effects in the current-current 
susceptibilities to see potential oscillations. 

Fermions in the confining model geometry clearly display bulk Friedel oscilla- 
tions. Because confinement leads to discrete bulk Fermi surfaces and each bulk Fermi 
surface is non-local in the radial direction, the bulk oscillations map into boundary 
Friedel oscillations. The loop parameter, "1/A^," sets the amplitude of the boundary 
Friedel oscillations, but not the wave vector at which they occur. 

For the electron star, on the other hand, the bulk Fermi surfaces are local features 
which vary along the radial direction. We find that the local character of the bulk 
Fermi surface prevents bulk Friedel oscillations from combining coherently. Each 
radial position sources a Friedel oscillation at a different boundary momentum. This 
is consistent with a boundary theory comprised of a continuum of Fermi surfaces, as 
proposed in [19]. 

Extrapolating from these two examples a stronger claim would be: For holo- 
graphic metallic states to exhibit a discrete set of boundary Fermi surfaces in all 
channels, there must be a discrete set of bulk Fermi surfaces and each bulk Fermi 
surfaces must be completely delocalized in the bulk. If either of these conditions fails 
any sharp features will be washed out. 

In this paper we only analyze systems with either completely delocalized bulk 
Fermi surfaces or completely localized surfaces. Testing whether the lessons learned 
from our two examples are general is of obvious importance. An alternative system 
is provided by the AdS black hole with Dirac hair [18] and it would be interesting 
to obtain the static susceptibility for this C3iSG clS well. It would also be interesting 
to repeat the analysis of fermions in a hard wall geometry with alternate boundary 
conditions. As discussed in [25], these should be dual to non- Fermi liquids and it 
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would be instructive to see how this simple change causes such drastic effects in the 
boundary theory. Finally, it would be instructive to incorporate the ideas of [25] in 
a more realistic model of confinement [33] and to see the effect of Fermi surfaces on 
other transport quantities. 
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A Perturbative Polarization Tensor 

In Section 2.1 we summarized the perturbative computation of the non-relativistic 
static density-density polarization, 

n(a; t, z' t') = -ie^G{x t, x' t')G{x' t\ xt) . (A.l) 

This is textbook material (see [31] for example) and details are only presented in this 
Appendix to make the paper more self-contained. 

As usual, it is simplest to analyze the polarization in momentum space, 

n(g, z/) = -le^ j j^^^g{k, uj)g{k + q,u + p). (A.2) 

Looking at the denominators of the time ordered Green's function (2.2) shows that 
only when k and A; + g lie on opposite sides oi kf will the answer be non-zero. 



n(g, ^) = e' j 



dk 



(27r)^ 
dk 



\q + k\- kf )d{kf - \k\) d{\k\ - kf)9{kf - |g + k\) 



V - {Ek+g - Ek)+ir] V - {Ek+g - Ek) - ir] 

-e{\q + k\-kf)e{kf-\k\) 



V - {Ek+g - Ek) +ir] u + {Ek+g - Ek) - if] J 



(A.3) 



In the second term we have changed momentum variables k —k — q and assumed 
that the spectrum is invariant under sign changes (that the Fermi surface is spherical). 
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The next step is to extract the imaginary part^ 

/dk 
—9(\q + k\-kj)9ikf-\k\) 

x[6{u- {Ek+q -Ek)) + 5{u+ {Ek+q - Ek))] . 

A.l Static Limit 

If we let — 7- with q fixed, we immediately see 

lmn(g,0) = 0. (A.5) 

There are two Fermi spheres in the problem. The step functions force you to be 
inside one and outside the other. On the other hand, the delta functions force the 
Fermi spheres to coincide and there is no solution to the constraints. 
From the static real part we have 

/dk 2 

Here we have used antisymmetry under k ■v^ k + q. In 2 + 1 dimensions we have 

ReU{q,0) = -^V [ dkdO eikf - k)—^— . (A.7) 
2vr^ J tjk+q - tjk 

For spherical Fermi surfaces the integrand is singular when q = 2kf. In general, the 
polarization tensor has non-analytic behavior in the momentum space and oscillatory 
behavior in position space. 

For example, for non-relativistic fermions we have Ek+q—Ek = q^ /2m+kq cos{6)/m. 
The polarization becomes, 

Ren(g,0) = / dkd9 9(kf - k)— ; — (A.8) 

' 2n^q J ^ ^ 'q/2 + kcos{e) ^ ' 

(A.9) 



2tx 




^Recall that for real 



^ ■p^^T7ri(5(a;-w')- (A.4) 



UJ ^ Uj' zk irj UJ — uj' 
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B Holographic Relative Polarization Tensor 

In Section 2.1 we introduced the relative polarization tensor, 

nZi^^^'^P) = - j ^^Tr {M^G{uXz,z')M^G{u,k + p,z',z)) 

-Tr (^M'^G°(u;, k, z, z')M''G\u, k + p, z', z)) . (B.l 
Using Eq. (3.4) we can perform the frequency integral: 



d'^k 



e{E,,{k + p))e{-E,{k)) e{-E,{k + p))e{E,{k)) 



Ee{k + p) - E,{k) + ir]s{E,{k)) Ee{k + p) - E,{k) - ir]s{Ee{k)) 

e{i')e{-e) e{-e)e{i) 



xTr XiA^)xi,iz')M'' Xi' ,k+piz')xl,^^{z) 

(Pk 
{2ny 



E'e,k+p - El, - iri E^,^^ - El - ir^ 

xTr xU^)x%{m'' • (B.2) 

In this we have denoted the contributions with zero chemical potential with ° and 
the primed summation indicates that the mode sums have no = terms. Also 
note that we have introduced the notation s{x) = sign(x). 

Though expression in Eq. (B.2) is quite complicated, there are two important 
features to note. Firstly, it is finite in the bulk UV. The large momentum and 
radial mode contributions do not know about the low energy finite density physics. 
Secondly, spherical bulk Fermi surfaces imply singular behavior due to the integrand's 
denominator aX p = 2kf and 1 = 1'. The only terms where 1 = 1' are when k and 
+ p lie on opposite sides of a Fermi surface. Quite generally, these singularities in 
the integrand will source Friedel oscillations in position space. 

B.l High Confinement Scale and Non-Relativistic Limit 

We will now consider a simplifying situation where there is a single Fermi surface 
occurring at /c = fc/ and ^ = 1. In this case we may split the summand in Eq. (B.2) 
into two classes, I = I' = 1 and all other combinations, 

(Pk ( e(\k + p\-kf)e{kf ~\k\) e{kf-\k+p\)e(\k\-kf) 



(27r)2 \Ei{k + p) - Ei{k) + ir]s{Ei{k)) E,{k + p) - Ei{k) ~ tr]s{Ei{k)) ^ 

xTr M^xi..(^)xI,,(^')M^Xi,fc+p(^')xU+p(^) + n^L,*.c- (B-S) 

If ^ = £' = 1 the integrand has a mild singularity due to the Fermi surface and there 
is no contribution from the vacuum terms. The rest of the mode number sums are 
analytic in the external momenta and do not involve the Fermi surface. 
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It is difficult to deal with the analytic contributions without explicit expressions 
for the energy eigenvalues and wave functions. However, we do note that their 
dominant contributions come from small momenta where the energy values scale 
with -^^(O) ~ Consequently, the analytic contributions to the relative vacuum 
polarization scale as Zm- If we model a high confinement scale, we can safely ignore 
all but the C. = C.' = 1 contributions, 

(fk ( e{\k + p\-kf)e{kf -\k\) e{kf-\k + p\)e{\k\-kf) 



(27r)2 \Ei{k + p) - Ei{k) + iris{Ei{k)) Ei{k + p) - Ei{k) - tr]s{E^{k)) , 

xTr XiA^)xUi^') Xi,.+p(^')xU+p(^) + 0{z^) . (B.4) 

While this is a great simplification, it is not quite enough for us to be able to 
find an analytic expression. In order to go further we will take a non-relativistic limit 
and scale the chemical potential towards the bulk mass. In this way we are left with 
a small kf and numerically observe that, for small momenta, the energy spectrum is 
well approximated by, 

Ei(fc) ~ ai + (B.5) 

Using this functional form the holographic calculation closely parallels the 2 + 1 
non-relativistic example in Section A.l with 2bi playing the role of an inverse mass. 



1 r (fk 



e{E,{k + p))e{-E,{k)) 



phi J {2TT)^lp + 2cos{e)k + ir]s{Ei{k)) 

xTr (m^ XiAz)xiki^')M'' xi ,fc+p(^')Xi,fc+p(^) 

e{E,{k-p))e{-E,{k)) 

p — 2 cos{9)k — iris{Ei{k)) 



Tr Xi,k-p{z)xlk-pi^')M^ XiA^')xlM 



(B.6) 

As in the non-relativistic case, the static polarization is purely real. We have 

iiZ{z,z',p} = -^v [ dkdo e{kf - k)——^ — — 

Ahin^p J p/2 + kcos{9) 

xTr (m^ XiMxiki^')M'' Xi,k+piz')xik+r>i^)) • (B.7) 

Numerically, it is expensive to work with Eq. (B.7) without further approximation. 
We can use the fact that kf is small which means that the magnitude of the loop 
momentum is never large and that the wave functions are slowly varying functions 
of the momenta. Therefore we may expand the wave function factors for small loop 
momenta, k. 



Tr (m° Xi,k{^)xU^')M'' XiMpi^')xU+,{z) 
~ Tr Xi,o(^)xI,o(^')M° XiA^IxIA^)) ■ (B.8) 
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This approximation loses some of the angular information in the wavef unctions, but 
is sufficient to display the essential features. In this way we find, 

nZ(z,z',p) = ^j^T^ {m" XiMxUz')M- xUz')xUz)) 

.vjikM0(k,-k)-^jL-^^. (B.9) 

Introducing the effective expansion parameter, A = -^^^ 

X - y^l-(^) % - 2fc/) j . (B.IO) 
This is the form of the polarization used in Section 3.4.1. 
C Flat Space QED Polarization Tensor 

We can understand the emergence of Friedel oscillations in relativistic Fermi liquid 
theory by studying how weakly coupled relativistic fermions respond to changes in 
an external gauge field. When working at finite chemical potential, the Lorentz 
invariance is broken and the polarization tensor may only be decomposed using 
the residual spatial rotation symmetry.^ If we introduce projections onto spatial 
directions, Pl, and the transverse direction, Pt, we may write 

W = GP^'' + FP^\ (C.l) 

Expressing the projectors in terms of the spatial momenta, we have 

P'^ = - (only spatial directions ^ 0) , (C.2) 

P^^ = Fr/A;^ - g^"" - P^^ . (C.3) 

As is mentioned in Section 3.3 in finite density computations, it is useful to 
further split the polarization tensor into a "vacuum" part and the relative, finite 
density, contribution, 

^"■^ = ^Zc + k:i , (C.4) 
iiZc = ^n,,T^o. (C.5) 



^For a general reference to relativistic field theories at finite temperature/density see 



-23- 



The vacuum contribution is a standard textbook computation [30] . Focusing on 3 + 1 
dimensions, F(0, k), may be found in [30, 37], 

'247r2\ „ , _ . .. o, / m 



e2 



(in \ 
— — 
fi + kfj 



k J \ y/W^^^^rri?kf - kfi 

A similar straightforward calculation yields 



e2 



4/i2 - 12m2 - gk'^\ f k - 2k f 



(4m2 + /c^)^/^ J y/W~T4m^kf - fik 



k \ a/P + Am'^kf + j2k 

The vacuum contributions contain the usual ultraviolet divergences for QED, and 
contain no information about the Fermi surface which is set by fi. For us, the most 
important feature is the logarithmic term with a singularity at = 2A;/ in both of 
these functions. Upon Fourier transformation this gives rise to spatial variations 
with wave number 2kf. When using this in the holographic calculation one should 
note that the appropriate values of k and kf are functions of the radial coordinate 
through redshift factors. 
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